On thermal fluctuations in quantum magnets 
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The effect of thermal fluctuations on the dynamics of a gapped quantum magnet is studied using 
inelastic neutron scattering on copper nitrate, a model material for the one-dimensional (ID) bond 
alternating Heisenberg chain, combined with theoretical and numerical analysis. We observe and 
interpret the thermally induced central peak due to intraband scattering as well as the thermal 
development of an asymmetric continuum of scattering. We relate this asymmetric line broadening 
to hard core constraints and quasi-particle interactions. Our findings are a counter example to 
recent assertions of universality of line broadening in ID systems and are to be considered as a new 
paradigm of behaviour, applicable to a broad range of quantum systems. 

PACS numbers: 75.10.Pq, 75.40Gb, 75.50Ee 



The behaviour of quantum systems at finite tempera- 
ture is of great importance for real applications as well 
as fundamental science. In spite of this few experimen- 
tal studies have hitherto been undertaken. The standard 
paradigm is thermally induced lifetime damping of the 
quasi-particles giving Lorentzian-type energy jTJ |5J [3] or 
wave vector jl] broadening. Indeed for one-dimensional 
systems linewidths have been proposed to display univer- 
sality [5j |6] . To test the applicability of these ideas more 
widely we studied a generic system of quasi-particles sub- 
ject to hard core interactions in the form of the gapped 
quantum magnet Cu(N03)2-2.5(D 2 0) by use of inelastic 
neutron scattering. We unambiguously identify different 
aspects of the interplay of quantum and thermal fluctua- 
tions. The most important finding is the thermal devel- 
opment of an asymmetric continuum of scattering, that 
differs strongly from the Lorentzian broadening familiar 
from conventional theories of thermal decoherence. We 
propose that this paradigm of behaviour is applicable to 
a broad range of low dimensional systems and that this 
key feature has been missed in previous work exhibiting 
the limitations on the parameter regime in which univer- 
sal behaviour may emerge. 

The alternating Heisenberg chain (AHC) of spin-1/2 
moments is described by the Hamiltonian H — JS^i • 
Si.2 + J'Si t 2 ■ <SWi,i where J and J' are antiferromagnetic 
exchange constants with J > J' and alternation parame- 
ter a defined as J' /J. The AHC has been an important 
paradigm in quantum magnetism for a long time [7J [5J [3] . 
Very recently it has attracted the attention of the quan- 
tum information community as an example for detecting 
and quantifying entanglement in solids |10j as well as en- 
coding and transporting q-bits |llj . 

In the AHC the dominant antiferromagnetic exchange 
interaction (J) between two neighbouring sites (sepa- 
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FIG. 1: a-b) One magnon excitations at 0.12K are resolu- 
tion limited. The dispersion and intensity are explained per- 
fectly by theory, (a) experiment, (b) theory. The wavevector 
modulation in intensity of the one-magnon band is from an 
interference term in the scattering amplitude from the two 
spins within the dimer. c) At 6K the one-magnon scattering 
has significant asymmetric broadening, and additional scat- 
tering appears around E=0meV. This intra-band scattering 
is found at low temperatures within the superimposed white 
boundary. 



rated by distance p and with chain repeat distance d) 
couples the spins into dimers, whose ground state is a 
singlet l/\/2{| tl> H !?>}• The dimer excited states 
are a spin-1 triplet with quantum numbers S z = 1,0,-1. 
To first approximation the global ground state of the 
chain is a product of dimer singlets. However, the intcr- 
dimer exchange (J 7 ) admixes a small amount (of order 
V3a/8) of polarised dimer pairs with total spin s = 
as ground state fluctuations. The elementary excita- 
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tions are a triplet of spin-1 states, or magnons. These 
are momentum k Bloch states obtained by exciting one 
dimer singlet to a triplet state. The inter-dimer coupling 
allows the excitation to hop from site-to-site along the 
chain. To lowest order in J' the one-magnon dispersion 
is uj(k) = J - (J' /2)cos(kd) [JJ]. 

Copper nitrate closely realises the AHC, the Cu 2+ ions 
have spin- 1/2 moments and the dominant exchange cou- 
plings are J=0.443±0.002 meV and a=0.227±0.005. For 
these parameters the magnon bandwidth is small com- 
pared to the gap. This parameter range is complemen- 
tary to the one analyzed in recent studies |4j [6] . Crucially 
the smallness of J'/ J leads to a clear separation in en- 
ergy between single magnon excitations and two magnon 
excitations even at high temperatures. 

We measured both the single and multiparticle states 
for a large (8g) single crystal [13] of 98% deuteration us- 
ing the OSIRIS neutron spectrometer at the ISIS Facility, 
UK |14j . At very low temperatures compared to the ex- 
change energy J, neutrons essentially measure the S 1 = 1 
spectrum by exciting from the ground state. The excita- 
tion to the one-magnon modes are easily seen and follow 
accurately the predicted cosinusoidal behaviour [1~4"1 H5]. 
The two-quasi-particle states can also be observed. These 
are sensitive to the interaction between particles and will 
be used later to quantitatively measure this. 

At low temperatures the dynamics of the strongly al- 
ternating Hciscnbcrg chain can concomitantly be thought 
of in terms of a low density gas of hard-core magnons. We 
follow the thermal evolution of correlations in this quasi- 
particle gas: The density of thermally excited magnons 
per-site can be estimated from the zero bandwidth limit 
to be n(T) w 3 exp(- J/kT)/{\ + 3 exp(- J/kT)). Mea- 
surements were made at temperatures T — 0.12,2,4,6 
K with densities of n = 6 • 10 -19 , 0.19, 0.45, 0.56 respec- 
tively (density at infinite temperature is n(oo) = 0.75). 
With increasing T the following picture emerges: the one- 
magnon mode develops into a scattering continuum, a 
new band of scattering emerges around zero energy with 
the same bandwidth as the one-magnon scattering, and 
the two-magnon scattering all but disappears, cf Fig. 1. 
Below we discuss each of these features in turn. 

(i) Temperature dependent lineshape of magnons. The 
effect of temperature on the scattering of magnons in 
copper nitrate was first investigated by Xu et al. |15j 
and discussed in the standard picture of a high tempera- 
ture (paramagnetic) phase, characterised by an absence 
of coherence, and a low temperature phase with quasi- 
particles, characterised by exponential decay, giving a 
Lorentzian lineshape. 

In sharp and remarkable contrast to this picture we 
find here magnon lineshapes which are highly asymmet- 
ric in energy. Whereas at T=0 the one magnon excita- 
tions are sharp modes that are resolution limited, com- 
parison to Figure 1 shows that the thermal evolution is 
characterised by the striking formation of a continuum 
of scattering weighted towards higher energies with in- 
creasing temperature. Experimental results are shown in 



Fig. 2 as cuts in energy at wavevectors ir/d and 2n/d. 

In a phenomenological description of these results we 
parametrise the lineshape by fitting it to a continuum 
with power law singularity at the lower (or upper) bound- 
ary uji (u> u ) convolved with a Lorentzian in energy of 
width r representing the damping, 



I(Q,w) 



;L/ r Q eQy-^(Q))eK(Q)-u/) (i) 



This is convolved with the instrumental resolution, mod- 
elled by a Gaussian with variance a = 0.012 meV. Re- 
sulting fits arc given as red lines in fig. 2(a,b), the widths 
of the Lorentzian and the power law exponent are shown 
in fig. 3(d,e). We argue below that strongly correlated 
behaviour from hard core constraints and magnon inter- 
actions explain these. 
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FIG. 2: Temperature dependence of intra-band and one- 
magnon excitations: The one-magnon scattering at a) 2n/d 
and b) ir/d, and c) the intra-band scattering at 57r/d is given 
for temperatures 120mK, 2K, 3K and 4K. The red dashed 
line is a parameterisation using equation 1, and the theory of 
interacting particles is shown as the solid blue line. Note that 
for 4K the agreement is less good because the series expansion 
is only valid at low particle density. 

(ii) Around zero energy a central peak in S(Q,u>) ap- 
pears due to intraband scattering within the one-magnon 
band: Neutrons can scatter via a change in momentum 
and energy (and possibly spin quantum number) off a 
thermally excited magnon. Precisely one scattering pro- 
cess exists for given wave vector Q and frequency w in a 
frequency range |a>| < ui m (Q) (where the maximum fre- 
quency u rn (Q = tt) is the bandwidth). To lowest order 
we have Lu rn (Q) = aJsin(Q/2) and from the density of 
states a square root at ±oj m (Q) results. In higher or- 
der this singularity becomes rounded out to a continuum 
[16] . Results for the lineshape of the central peak are 
displayed in Fig. 4(c). Cuts through the intraband scat- 
tering for Q = 7r for different temperatures are shown in 
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Figure 3(c). Although the central regime is masked by 
other incoherent scattering processes from the material, 
the increase towards u m and the drastic decrease beyond 
uj m is clearly seen. 




b) 


1 Q=7t/d 




J{ T=2K 




J\ T=3K 




/\ T=4K 





0.4 0.5 0.5 

Energy (meV) 




u " #. e) 



T(K) 



FIG. 3: Asymmetric lineshapes as a function of temperature. 
The theoretical lineshapes of the one-magnon band at a) 2%/d 
and b) 7r/d, and of the intra-band at c) 5n/d, calculated using 
the theory of interactiing particles described in the text, d) 
and e) give the fitted linewidths and exponents respectively 
as described by equation (1). 

The basic square root singularity characterising this 
scattering process appears similar to the scattering from 
excited domain walls in quantum Hciscnbcrg-Ising chains 
as predicted by Villain [17] and observed later in the 
quasi ID magnets CsCoCl 3 [HI US] and CsCoBr 3 [20] . 
However, this similarity is to large extent formal (cosine 
spectrum of the basic excitation) whereas the basic pro- 
cesses behind the scattering are very different [21] : In the 
AHC the leading contribution to the intraband scatter- 
ing at low temperatures arises from scattering processes 
between two well defined single particle states. On the 
other hand, spin excitations in the Heisenberg-Ising chain 
occur in pairs of domain walls and the leading contribu- 
tion to the low temperature intraband scattering is due 
to scattering on one member of the domain wall pair. 
Whilst the topological nature of these soliton-like excita- 
tions restricts it to ID, the dimer case studied here can 
serve as a paradigm for gapped quantum magnets in any 
number of dimensions. 

(iii) Two-magnon scattering and bound states: Neu- 
trons scatter from the weakly admixed polarised dimer 
pairs composing the ground state fluctuations to the two- 
magnon states. This scattering is far more sensitive to 
the composition, interactions and phase factors between 
particles than the one-magnon signal and so provides key 
insight. In particular we resolve for the first time the 
spin-1 bound mode of the AHC which is a direct signa- 
ture of the two particle interaction, see Figure 4. 

The excitations obey a hard core constraint as each 
dimer site can be occupied by at most one magnon 
at a time. In addition to this particles occupying 
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FIG. 4: The two-magnon continuum and bound mode at 
T=0.12K. a) the superimposed line gives the upper and lower 
boundaries of the two magnon scattering and the position of 
the bound mode, lying just below the node, b) is a cut of the 
scattering intensity at Q = 2n/d as a function of energy, scat- 
tering is observed between the upper and lower two-magnon 
envelope, c) is a cut at the node position (Q = 37r/d). The 
node occurs at an energy of 0.88meV, however it is predicted 
to have very little intensity at this wavevector. Red lines are 
calculations based on the two particle Schrodinger solution 
[22] . A bound mode is observed at the theoretically predicted 
lower energy of 0.855meV. 



neighbouring sites interact with each other through the 
inter-dimer exchange (J'). This interaction only exists 
for nearest neighbours and is absent for larger dimer 
separations. The nearest neighbour potential, V Sl is 
— J'/2, — J'/4, J'/4 for pairs with combined total spin 
s = 0,1,2 respectively and provides an attraction that 
has important effects on multi-particle states. There 
are two types of solution to the two-particle Schrodinger 
equation [55]: The first has the form of a continuum of 
states, where fci and k 2 are momenta of the magnons 
along the chain, giving a total momentum Q — k\ + k 2 
and energy tkj — hu>(k 1 ) + haj(k 2 ). The scattering ma- 
trix for two particles in a state with total spin s is given 

by S s kiM = ~ } + Cul]ttl\*-^C^ , where A n = 0, 
Ai = 1/2, 
shift is 



l+ e -i(ki+k 2 )d_ 2A s e~ ik i d ' 

A2 = —1/2, and the corresponding phase 
b%^ k2 = ~l\n{S s ki M ). We note that in general 
the S-matrix differs from the non-interacting hard core 
value Sk lt k 2 = — 1 but reduces to the latter in the limit 
of vanishing momenta. 

The second solution to the two-particle Schrodinger 
equation has the form of an exponentially decaying 
bound mode which exists over all Q for s — and around 
the nodes in the continuum (Q ~ ft/d) for s = 1. For 
s = 2 the mode occurs at energies above the two-magnon 
continuum ( "anti-bound mode" ) . The binding energy 
and intensity of the bound mode for low temperature 
agree with the calculations [25] as do the deconfined two 
particle states, see Figure 4. This experimental observa- 
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tion of the bound mode reveals unambiguously the pres- 
ence of interparticle interactions. 

The decrease in the two-particle scattering with in- 
creasing temperature is a direct result of the hard- 
core constraint. A neutron cannot create a further two 
magnons on a given pair of dimers if any of the dimcrs 
is already occupied by a magnon. The density of adja- 
cent dimer pairs that are both unoccupied falls approxi- 
mately as n ss w l/(l + 6exp(-J/fcT)+9exp(-2J/fcT)), 
i.e. n ss = 1,0.47,0.2,0.14 respectively for temperatures 
0.12, 2,4, 6K and reaches a value of n ss (oo) — 0.0625 at 
infinite temperature. The intensity of the band is 1 % of 
the total scattering at lowest temperatures [22] and with 
an intensity of only 0.06 % at high temperature becomes 
very hard to detect. 

To understand further the asymmetric thermal line- 
shapes for the magnons we analyzed the T > dynami- 
cal response in the framework of a low density expansion, 
where the small parameter is e~ J l T |16j . In this regime, 
lineshapes agree with those calculated by full diagonali- 
sation of finite chains [23]. From these diagonalisations 
one also concludes that the asymmetry persists beyond 
the low temperature regime of the analytical approach, 
at least up to temperatures of twice the gap. 

The asymmetry of the lineshapes results from the com- 
bined effect of the hard core and additional interactions 
between neighbouring sites, and the density of states: 
For hard core particles with no contact potential, i.e. 
V = 0, the quasiparticlcs are fcrmionic in the sense 
that their scattering matrix is -1. The broadening of the 
magnon line is found to be highly asymmetric in this case. 
The presence of the contact potential (" stickiness" ) be- 
tween quasiparticles further modifies the lineshape. This 
demonstrates that, as expected on general grounds, the 
broadening of the line reflects the details of the quasipar- 
ticle interactions. On a qualitative level the broadening 
of the magnon line can be understood by considering the 
joint density of states for transitions between thermally 
occupied one magnon states and unoccupied two magnon 



states 

P,Pl,P2 

Here n(jp) is the thermal occupation number for a one- 
magnon state with momentum p and n(pi,p2) is the 
probability that the two-magnon state characterized by 
momenta pi,p2 is unoccupied. For —Tr/2d < Q < n/2d 
and it i '2d < Q < "in/ 2d this function is skewed towards 
higher and lower energies respectively, in agreement with 
experimental (Fig. 2 a,b) and numerical lineshapes. The 
specific form of the lineshape results from the matrix ele- 
ment and hence by the magnon-magnon interaction, the 
lines in Figures are such calculations. 

This behaviour then represents a paradigm for ther- 
mal effects in quantum magnets which differs sharply to 
the picture of simple lifetime broadening due to thermal 
decoherence as is observed for example in the 2D Hciscn- 
berg antiferromagnet [3]. We expect that such strongly 
correlated physics is generic to ensembles of interacting 
hard core particles and should be seen in other quantum 
gases including other gapped quantum magnets, and in 
higher dimensions. 

In summary, we have performed a comprehensive study 
of thermal effects in a quantum magnet. Our study shows 
that temperature does not give rise to simple decoherence 
but instead promotes the formation of a strongly corre- 
lated gas of quasi-partcles. A description of this state has 
been achieved using non perturbative methods. The find- 
ings here present a model of such behavior that should 
apply to a large range of quantum systems. 
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